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Abstract 

The present observations of Cabbibo angle and solar mixing angle satisfy the empirical relation 
called Quark-Lepton Complementarity (QLC), namely 9[ 2 ~ vr/4 — 9q. It suggests the existence of 
correlation between quarks and leptons which is supported by the idea of grand unification. We 
propose a specific ansatz for the structure of Yukawa matrices in SO (10) unified theory which 
leads to such relation if neutrinos get masses through type-II seesaw mechanism. Viability of this 
ansatz is discussed through detailed analysis of fermion masses and mixing angles all of which can 
be explained in a model which uses three Higgs fields transforming as 10 and one transforming as 
126 representations under 50(10). It is shown that the proposed ansatz can be derived from an 
extended model based on the two pairs of 16-dimensional vector-like fermions and an £4 flavor 
symmetry. The model leads to the lepton mixing matrix that is dominantly bimaximal with 
O(0c) corrections related to quark mixings. A generic prediction of the model is the reactor angle 
d[ 3 ~ 6c /y/2 which is close to its present experimental upper bound. 
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I. INTRODUCTION 



Experiments on neutrino oscillations have revealed that two of the three leptonic mixing 
angles are large. One of them called the atmospheric mixing angle is almost maximal 
#23 = 42.3° (i^VJ and the other known as the solar mixing angle 9[ 2 = 34.5° (jj s) i s smaller 
compared to it I n contrast, the observed quark mixing angles are small and hierarchical. 
The largest angle is the Cabibbo angle 9c = 9 US « 13° while other two are 9 c b ~ 2.4° 
and 8 u b ~ 0.2°. An understanding of such wide dissimilarity between the quark and lepton 



mixing patterns is considered as one of the ma 



standard model. It has been observed long ago [2] that there exists an interesting empirical 
relation between quark and lepton mixing angles. 



or challenges for the physics beyond the 



0[ 2 + 0us~l (1) 



The above relation is known as Quark-Lepton Complementarity [3|-|6[ and still favored by 
the present experimental data within their measurement errors. It is also possible to write 
similar relation between 23 angles of quark and lepton mixing. 

9 l 23 + e cb ~^ (2) 

If such relations are not accidental, they strongly suggest the common roots between quarks 
and leptons jsrjsj. Clearly it is very hard to realize such relations in ordinary bottom- up 
approaches where the quarks and leptons are treated separately with no specific connections 
between them. So one requires top-down approaches like the Grand Unified Theories (GUT) 
which sometime also unify quarks and leptons and provide a framework to construct a model 
in which QLC relation can be embedded in a natural way. 

The general conditions under which QLC relation ([1]) can be realized from quark-lepton 
unification are thoroughly discussed in j^, 0]. We describe one such possibility here. The 
quark mixing matrix known as Cabibbo-Kobayashi-Maskawa(CKM) matrix is defined 
as Vqkm = U}jJd where U u (Ud) is unitary matrix that diagonalizes the up- (down-) type 
quark mass matrix. Corresponding leptonic mixing matrix, also called Pontecorvo-Maki- 
Nakagawa-Sakata(PMNS) matrix, is Vpmns = U\\J V . Assume that the structure of neutrino 
and quark mass matrices at high scale are such that the PMNS matrix is exact bimaximal 
Vpmns — Ubm whereas the CKM matrix is an identity matrix to a leading order. 
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Both the mixing matrices get corrected by O(0c) terms coming from the next leading order 
where the down quark and charged lepton mass matrices are equal (or nearly so) . In this sce- 
nario, a QLC relation can emerge from quark-lepton unification at high scale. Construction 
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of a realistic GUT model in which all fermion masses and mixing angles are correctly repro- 
duced along with QLC is highly non-trivial. In fact several models |5| proposed to explain 
QLC are based on a smaller gauge group, namely Pati-Salam S77(4) c x SU{2) L x SU(2) R 
group. A complete and realistic model based on 50(10) GUT has not been proposed so far. 
The original proposal 0] was based on SU(5) relation M e = Mj but detailed explanation 
of the fermionic spectrum was not developed. Here we present a predictive 50(10) based 
unified description of fermion masses and mixing in which QLC relation can be naturally 
realized. 



The renormalizable models based on the 5*0(10) gauge group are quite powerful in 
constraining the fermion mass structure. Moreover, they provide a natural framework for 
understanding neutrino masses because of the seesaw mechanisms inherent in them. Fermion 
masses arise in these models through their couplings to Higgs fields transforming as 10, 
126 and 120 dimensional representations under SO (10). Neutrino masses arise from two 
separate sources either from the vacuum expectation value (vev) of the right handed triplet 
(type-I) or from the left handed triplet (type-II) Higgs. The minimal model with 10 and 
126 Higgs fields has attracted a lot attention [7H9]. There also exist a class of models where 
appropriate flavor symmetry is integrated with 5*0(10) framework with extended Higgs sec- 
tor 10Nl3l] to construct a predictive theory which can simultaneously explain hierarchical 



nature of quark masses and mixing angles and large lepton mixing angles. In this paper, we 
show that QLC relation follows in a specific 5O(10) model combined with 54 symmetry if 
dominant source of neutrino mass is type-II. An additional Z n symmetry is required in the 
model to get desired interactions between various fields. 



The paper is organized as follows. We describe the fermion mass relations in the model 
based on renormalizable supersymmetric (SUSY) 5O(10) GUT in the next section. In sec- 
tion III, we propose a specific ansatz which predictively interrelates various observables of 
quark and lepton sectors and leads to QLC relation. We also discuss the phenomenological 
implications of such ansatz in this section. In section IV, we justify the proposed ansatz by 
a flavor symmetry group S4 x Z n . Finally we conclude in section V. 



II. RENORMALIZABLE SUSY 5O(10) MODEL FOR FERMION MASSES 



We consider three families of 16-dimensional fermions obtaining their masses from renor- 
malizable couplings to four Higgs multiplates, three of them (denoted by $, $' and $") 
transforming as 10 and the other (£) as 126 dimensional representations under 5O(10). 
The 5O(10) breaking can be achieved with 210 + 54 + 126 + 126 [9]. The Yukawa inter- 
actions of the model can be written as 



10 



(4) 
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where Y> are symmetric Yukawa coupling matrices. The representations $, <&" and S 
have two minimal supersymmetric standard model (MSSM) doublets in each of them. It is 
assumed that only one linear combination of the up-type doublets and one of the down-type 
doublets remain light and play the role of H u and H d fields. Once these light doublets 
acquire vacuum expectation values, they break electroweak symmetry and generate the 
fermion masses as well. The resulting fermion mass matrices can be suitably written as 

M d = H + F + t H' + H"; 
M u = rH + sF + H'+p H"; 
M e = H~3F + t H' + H"; 
M D = r H-3s F + H'+p H"; 
M L = r L F; 

Mr = r- R l F (5) 

where H, F, H' and H" are obtained by multiplying electroweak vevs and Higgs mixing 
parameters with Yukawa coupling matrices Y w , Y 10 > and Y w n respectively. r,s,t,p,r L 
are tr are dimensionless parameters determined by the Clebsch-Gordan coefficients, ratios 



of vevs, and mixing among Higgs fields (see 11| for example). denotes neutrino Dirac 



mass matrix. Ml(Mr) is the Majorana mass matrix for left-(right-)handed neutrinos which 
receives a contribution only from the vev of £ field. In generic SO(10) models of this type, 
the effective neutrino mass matrix Ai y for the three light neutrinos has type-I and type-II 
contributions. 

M v = M'J + Ml = r L F - r R M D F- x M T D . (6) 

In general, both contributions are present and they depend on two different parameters so 
one may dominate over the other. It has been shown in several references [o| that it is 
possible to have symmetry breaking pattern in SO (10) where type-II term dominates over 
the type-I contributions. In this limit, neutrino masses and mixing are governed by F which 
can be written as F ~ M d — M e . It is well known that this relation establish interesting 
relationship between b — r unification and large atmospheric mixing angle The equations 
([5]) and ([6]) are the key equations that provide basic platform to construct a model in which 
the QLC relation (pQ) can be realized. 



III. ANSATZ 

We propose following ansatz which leads to relation ([!]). 
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where I is 3 x 3 identity matrix. To do the simple analytical study of such ansatz we assume 
that all the above parameters are real. Without loss of generality, we can express the above 
matrices in a basis with diagonal H. Such basis are obtained by rotating the 16- dimensional 
fermion fields in 2-3 plane by an angle it/ 4. The matrices in ([7|) will be redefined in new 
basis as 

(H, F, H', H") R 23 (£) (H, F, H>, H")R T 23 (£) (8) 
and can be rewritten as 
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Before we present the detailed analysis let us look at some immediate implications of the 
above ansatz. The dominant 10-Higgs coupling matrix H has rank-1. As it was pointed out 
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this can simultaneously explain both the observed hierarchy of quark masses as 
well as the origin of large lepton mixings if the light neutrino masses are generated through 
type-II seesaw mechanism. Assuming only one 10-Higgs H contribution in charged fermion 
mass matrices, we get at zeroth order, 
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(10) 



Correct b — r unification and large lepton mixings (bimaximal) are obtained with no mixings 
between quarks. The charged fermions of first two generations are massless in this case. 
Further, the contributions coming from other Higgs coupling matrices F, H' and H" make 
the model realistic by giving nonzero masses to first two fermion generations as well as by 
perturbing both the mixing matrices which reproduce observed mixing patterns for both 
the quark and lepton sectors. 



We now present the detailed analysis of ansatzQ. Substituting it in eq.flS]) and (jSJ), we 



get 



' s(b + c) + x 1 sa 



a 



sa + a' \ 
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sa — a' sb + x' sc 
sa + a' sc rh + sb + x' J 

^ —3(6 + c) + x —3a — ta! —3a + ta' ' 

—3a — ta' —3b + x —3c 
^ —3a + ta' —3c h — 3b + x j 



Md 



^b + c + xa — ta' a + ta' ^ 

a — ta' b + x c 
^ a + ta' c h + b + x j 



r L F 



(11) 



where x' = px. Since each mass matrix is real symmetric, it can be diagonalized by a 
rotation matrix parameterized (in the standard parameterization) by three angles. 



R T f M f R f = Diag(m fl ,m f2 ,m f3 ); 

Rf = R2^{ 3 )Rl2,{0[z)Rl2{d[ 2 ) 



(12) 
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where / = d,u,e,v and i?y is a rotation matrix in ij plane. The charged fermion mass 
matrices are hierarchical (h 3> b, c 3> a, a' ^> x, x') and can be approximately diagonalized 
by Jacobi rotation. The results obtained from such diagonalization for the quark sector are 
displayed below. 



( c 2 
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m ^ fe + c + x _ v (13) 



m t « r/i + sb + x' + (9 
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+£i£ v 9 __c ; 9 , /3S3 _^: (l+e . (15) 



sa — a' /„ sc — x'\ „„. sc _ sa + a' 



( 2 + — j ; l 1 + wJ' (16 » 

Let us underline some important points in connection with above relations. 

• The six real parameters h, 6, x, r, s, x' can be approximated from the six quark masses. 
rrif, and m s determine the parameters h and 6. It is easy to see that r « mt/mi, and 
s « m c /m s are required to obtain the masses of heavy quarks m ( and m c . Further, 

and m u fix the values of x and x' . Since 6, c ^> x, we require c ~ — b to obtain 
small masses of first generation fermions. 

• Let us assume that a' ~ sa in order to keep ^ ^12 • Also note that #23 ~ (s/r)8 d 3 <C 
6*23 and 0" 3 ~ ( s / r )^i3 *C #f 3 . In this limit, the quark mixing matrix takes the form 

Vckm = UtU d &U d & R 23 (d d 3 )R 13 (d d 13 )R 12 (8 d 2 ) (17) 

• The elements of the CKM matrix fix some more parameters as follows. 

c ~ — V c b h] a — ta' ~ — V us 6; a + ta ~ — V u b h. (18) 
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An interesting relationship between V us and V u b can be found in the limit t ~ 0. 

Kb »K,- + f—fl (19) 

We will show later in this section that t ~ is a necessary requirement to obtain QLC 
relation([T|). 

• Our assumption of real parameters makes the theory CP invariant. The observed 
CP violation in the quark sector can be accommodated by making some parameters 
complex. 



It is interesting to note that all the parameters are fixed in terms of the observables of 
the quark sector. Hence the entire lepton sector emerges as the prediction of the model. Let 
us first derive the predictions for the charged leptons. 

m T « h - 3b + x + O ( — 

(3a + ta') 2 / x\ . n fc 2 

(3a + ta') 2 / x\ _ f a 2 \ 
m .«-S(l +e ) + x + ^- 5r L(l + _) + 0^j. (20) 

3a + ta' ( 3c+x\ „ 3c 3a — ta' / c\ . . 

( 2 + — > — + (21 » 

Noteworthy features of the above relations are the following, 

• It predicts m T m m& and m M w — 3m s . 

• For 6 = — c, m e « which is viable with observed values of m e and extrapolated 



at the GUT scale within 3a deviations 15] . However for b ^ — c, any desired value of 
md/m e can be obtained. 

For t ~ 0, 0f o « 0^, 0?, ~ -30,* and 0?, « -39, lh . 



The light neutrino mass matrix in eq. (fTT|) has the most general form which can be diag- 
onalized by bimaximal matrix Ubm- The mass eigenvalues are, 

mi = mo(b + c + \/2a); m 2 = m (b + c — V2a); m 3 = m (b — c) (22) 

Interestingly, for b = —c (which can now also be written as ~ m s /m(,), we get the partial 
degenerate neutrino mass spectrum ni\ = — m2 <C which leads to vanishing solar (mass) 2 
difference (Am 2 oi = m 2 — m 2 = 0) at high scale. We performed numerical study and found 
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that the radiative corrections to the original neutrino mass matrix are unable to generate 
the required splitting between mi and m^- Another way to induce non zero value of Am 2 sol 
is to allow type-I contribution to the original type-II seesaw neutrino mass matrix. However 
such contribution is highly hierarchical (like Mf) and it largely contributes to the 33 element 
of neutrino mass matrix which ultimately spoils the nice symmetry of neutrino mass matrix 
and hence the bimaximality of neutrino mixings. This forces us to consider the case where 
Vcb 7^ m s / m b- In this case we obtain the following expression for the ratio of the solar to 
atmospheric squared mass difference. 

AmL ~ V2V US ( ^ + (23) 

atm\ r ^ J 



Am 2 atm V v cb V m b 



Note that one requires m s /m h ~ 1.08 V c b to obtain the observed value of Am^/Am^ 
0.031) and it implies that V c b < m s /mb which is not favored by their present observed values 



extrapolated at the GUT scale. However as argued in [13j, the threshold corrections to b — s 
quark mass mixing from gluino and wino exchange via one-loop diagrams can give desired 
value of V c b- The required deviation from b = — c is quantified by 

b + c ^ m s ( 1 — ) < 0.08 m s 

\ m s /m b J 

which is small and of order of first generation fermion masses and hence allows the correct 
ma in eq.(|T3l). 



The leptonic mixing matrix can be seen as dominant bimaximal mixing resulting from 
neutrino mass matrix and then corrected by 0(9q) terms coming from the unitary matrix 
U e which diagonalize charged lepton mass matrix. 

V PMNS = U\U V = UjU BM (24) 

where U e = i?23(— 36 l c fe)-Ri 3 (— 39 ub )Ri 2 (9c)- The resulting neutrino mixing parameters are 
the following. 



TJ -(V ) ~ 1 + ^- 3 ^) 

Ue2 — \ VPMNS)12 r 7= T 



V2 2 

= (VpmNs)23 » + 3 ^&); 

U e3 = (V PMNS ) 13 « -4|(K S + 3K b ). (25) 

The correction of 0($c) from charged lepton generates correct solar mixing angle 
which follows QLC relation (CQ). The atmospheric mixing angle gets considerable devia- 

7T 

tion #23 ~ ^ + 39 c b in this model unlike the standard QLC relation for 23 mixing angle of 
quark and lepton given in eq.(T5]). The model also predicts large value of U e -s ~ 0.16 which 
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can be tested in planned long baseline experiments. 



Note that eq. fl25l) holds at GUT scale which might be changed by RGE corrections in 
principle. However it is known that running of the Cabibbo angle is negligibly small in 
MSSM even with large value of tan/3. Running of leptonic mixing angle depends on the type 
of mass spectrum of light neutrinos. For b ^ — c, neutrino mass spectrum follows normal 
hierarchy m x < m 2 <C m 3 . The effect of RGE corrections are known to be negligible in this 
case and eq. (l25l) holds at low scale also. 

We now provide an example of values of the parameters of eq.(JH]) which successfully 
generate entire fermion mass spectrum as well as mixing patterns for both quark and lepton 
sector. The required CP violation in the quark sector is incorporated by making a' complex. 
In the limit t ~ 0, a' contributes only to the up quark mass matrix and does not change 
the other predictions of ansatz given in eq-Q- One more parameter x' is made complex to 
reproduce m u correctly. The numerical values of parameters are 

h = 1.7 GeV; b = 0.0243 GeV; c = -0.022113 GeV; a = -0.0052 GeV; 
a' = (0.0344247 - 0.028885z) GeV; x' = (0.0233596 - 0.00293374i) GeV; 
x = 0.00325 GeV; r = 55.88; s = -8.64198; t = 0. (26) 

Substituting these numbers in eq. ffTll . we get 

m t = 94.8 GeV 
m b = 1.73 GeV 
m T = 1.63 GeV 



m c = 0.19 GeV; m u = 0.65 MeV; 
m s = 28.5 MeV; m d = 4.21 MeV; (27) 
m. = 75.4 MeV; m e = 0.35 MeV. 



sm6 us = 0.222; sinfl c6 = 0.015; sm6 ub = 0.005; 8 C km = 60.9°; 



(28) 



sm 2 9[ 2 = 0.368; sin^< 3 = 0.527; sin% = 0.024; = 0.030. 

The obtained spectrum is in good agreement with the data extrapolated at the GUT 
scale. For example, we compare our results with the charged fermion masses obtained at 
the GUT scale in the MSSM for tan/3=55, M SUS y = 1 TeV and M GUT = 2 x 10 16 GeV given 
in table 5 of reference Q. All charged fermion masses (except md) obtained here fits with 
the data within la. Our ansatz predicts larger value of m^. The quark mixing angles 9 cb 
is small (< m s /m&) as required by eq.(l23l. The reproduced values of lepton mixing angles 
and Aml ol / Am atm are also in accordance with their updated low energy values (within 3cr 
measurement errors) given in [l|. 



IV. THE MODEL 

In this section, we will illustrate how the ansatz (CO) can be obtained in a model from 
flavor symmetry. We use discrete flavor symmetry based on the group 64 which is a group 
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of permutation of four distinct objects. It has 24 distinct elements filled in five conjugate 
classes and hence five irreducible representations of dimensions 3 2 , 3i, 2, 1 2 and li. A sin- 
glet representation with subscript "2" changes sign under transformation involving the odd 
number of permutations of S4. More details on the grou p th eory of S4, its multiplication 
rules and the Clebsch-Gordan coefficients are reported in 14 



Our model follows the same line as model constructed in [13] and uses the same sym- 
metry group. However it differs at some places since the ansatz required here is different 
from their ansatz. The basic matter fields and Higgs fields content of the model is the same 
as discussed in section II. In addition to this we use five flavon fields which are singlets 
under SO (10) and two pair of vector-like fermion fields which transform like 16 © 16 under 
50(10). We impose the S4 symmetry together with Z n symmetry to get desired structure 
of Yukawa matrices. Three matter fields ip are assigned as 3 2 dimensional representation of 
5 4 while five flavon fields x^,r\^a and a' form 3i,3 2 ,3i,li and 1 2 representations of 5 4 
respectively. The other fields are singlet (li or 1 2 ) under 5 4 . An additional Z n symmetry 
is required to allow/forbid interactions between particular fields. The Z n charges of various 
fields are listed in tabledU) where ui = e^ 27T ^ n \ 
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TABLE I. Various fields and their representations under 5O(10) x 5 4 x Z, 



Let us consider a theory above GUT scale which is invariant under the symmetry group 
5O(10) x 54 x Z n . The Yukawa superpotential allowed by such symmetry can be written as 



W = (4>1p)^Vl + X^Vl^Vl® + Mi^yi^yi 

+ E ^(X 2 ^),S + + -^<x 2 (^)X 

i 

a' a" 
+ -a'( v ^' + ^a' 2 (^" (29) 

where A is the Planck scale up to which the theory is valid. The 5 4 singlet contraction of 
flavor index is indicated with bracket, ccj, a', a", (3, 7, A, and A' are coefficients of 0(1). The 
term (x 2 ipip) represents all the different 5 4 contractions which can be constructed as follows: 



(X 2 ^)i = ((XX) ii(#)ii), (to)a(#)2), ((xx)si W)ai), ((xx) a a W)s a ), 

((x^MxVOO, ((xV>)2(xV>) 2 ), ((x^MxVOO, ((x^Mx^aJ (30) 
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where (...)# indicates the representation R under S4. Now consider a theory below the scale 
of Mi t 2 and at the GUT scale. The effective superpotential after integrating out heavy 
vector-like fields is given by, 

w eff = -^{u){H)<b + ILwmje 



Ml 
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(31) 



where first two terms allow the desired rank-1 structure of Yukawa matrices. Note that 
effective Yukawa superpotential still has the symmetry SO (10) x S4 x Z n . This symmetry 
will be broken to SO (10) by vevs of the flavon fields. In order to get the desired structure of 
Yukawa couplings, we will choose particular vacuum alignment of the flavon fields as given 
below. 
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These vevs of flavon fields break flavor symmetry S4 at the GUT scale and generate following 
structure of various Yukawa couplings. 
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where all non-relevant Clebsch-Gordan coefficients are suitably absorbed, a and d Q are 
linear combinations of different c^. The Yukawa matrices derived from the super poten- 
tial can successfully explain the ansatz given in eq.((7]). Note that M2 <C A which implies 
b + c — c (or 6 ^ —c) in eq.([7]). Further, the assumption Mi M 2 leads to h 3> b, c. 
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It is very important to show that the required vacuum structure of flavon fields (132]) is 
allowed by flavon superpotential. This point has already been discussed in great details in 
reference 
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Since our model has the same kind of flavon structure as theirs, we simply 
use their results. Note that due to non-trivial Z n charges, bilinear terms which correspond 
to masses of flavon fields are not allowed. As a result of this the model requires doubling 
of flavon fields to allow Dirac type mass terms. The new flavon fields have the same S4 
representations but opposite Z n charges. It has been shown in [13j that all the desired vacua 
of eq.( )32|) are present in the model. 



V. SUMMARY 

In this paper, we have studied a possible way to realize QLC relation ([T|) between the 
Cabibbo angle and solar mixing angle in realistic quark-lepton unification theory based on 
SO (10) gauge group. We have shown here that it is indeed possible to obtain such relation 
starting from the fermionic mass structure §5§ if they are supplemented with ansatz ([7]) and 
assuming that only type-II seesaw mechanism is responsible for light neutrino masses. One 
necessary ingredient for QLC is bimaximal mixing pattern from the neutrino sector which 
has been obtained through specific ansatz. Our ansatz also makes use of recently proposed 
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rank-1 strategy which naturally explains charged fermions mass hierarchy as well as 
origin of hierarchical quark mixing angles as opposed to the large lepton mixing angles. 
We have shown through the detailed analysis that this ansatz is capable of explaining the 
entire fermionic spectrum and not just the QLC relation. Moreover, the various predictions 
made by such ansatz are in agreement with observations. We have shown that the proposed 
ansatz can be obtained in a model from a discrete flavor symmetry group S4 together with an 
additional Z n symmetry. A generic prediction of our approach is 9[ 3 ~ 9 C / a/2 which is near 
to its current experimental upper bound. The atmospheric mixing angle gets considerable 
deviation from maximality (#23 ~ 7r/4 + 30 c b) in this approach. These predictions can be 
confirmed or excluded by the current generation of neutrino oscillations experiments. 
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